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Abstract. In this paper, we consider the following nonlinear eigenvalue prob- 
lem for the Monge-Ampere equation: find a non-negative weakly convex clas- 
sical solution / satisfying 

!detD 2 f = fP inQ. 
f = tp on dQ. 

for a strictly convex smooth domain Q C R 2 and < p < 2. When {/ = 0} 
contains a convex domain, we find a classical solution which is smooth on 



{/ > 0} and whose free boundary d{f = 0} is also smooth. 



1. Introduction 



We consider in this work the following nonlinear eigenvalue problem for the 
Monge-Ampere equation: find a non-negative weakly convex classical solution / 
satisfying 



(MA) 



detD 2 f = fP in ft 
f = (f on 9f2 



for a strictly convex bounded smooth domain fi C R 2 , with cp > on <9f2 and 
smooth, and < p < 2. 

The study of problem (MA) is motivated by the general Minkowski problem in 
differential geometry, asking to find the manifold whose Gauss curvature has been 
prescribed. More generally, the Gauss curvature itself may depend on the graph 
z = f(P) of the manifold, namely 

det D 2 f(P) = h(P,f(P)yf(P)). 

For a positive bounded h, this problem has been discussed by many authors and 
the C 1, ^regularity of / has been established (c.f. in |GTj ). When h is allowed to 
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be zero, / is not always a C 1 ' function, as it will be discussed in the sequel. The 
regularity of / is an open problem (c.f. Aubin j Au] ) . 

One of the interesting cases is when ft = on the vanishing set of /, especially 
when h = f p , as in problem (|MAp . For p = 1, this problem corresponds to an 
eigenvalue problem describing the asymptotic behavior, as t — > T, of the parabolic 
Monge- Ampere equation 

f t = detD 2 f inOx(0,T) 
/ = rj(t) <p on dVl 

where t]{t) = l/(T-t). 

For / < 0, / = on d£l and h = (-/)" (n+2) in Q, problem (|MA"f was considered 
by Cheng and Yau in CY . When, h = (—A/)™ problem (|MA[1 corresponds to 
the eigenvalue problem for the concave operator (det D 2 f)i and has been studied 
in [Li] . The exponential nonlinearity, h = e -2 ^ has been studied by Cheng and 
Yau in |CY2j . Equation detD 2 / = h with a degenerate source term h has been 
studied at [G]. The limiting case p — > + , f(x) p — > X{/>o} was considered by O. 
Savin, [S] as the obstacle problem for Monge- Ampere equation, where the obstacle 
stays below the graph of /. The second author also considered the case where the 
obstacle stays above the graph of /. 

Since f(x) p — > X{f>o} as p — > + , (IMAI) corresponds to a perturbation problem 
for the obstacle problem 

det D 2 f = X{f>o} 

and / > in ft. 

Depending on the boundary values if one of the three possibilities may occur in 
(MA): 

i. / > in f2: the equation (MA) is then strictly elliptic and by the regularity 
theory of fully- nonlinear equations, / is C°° smooth in il (cf. |CCj ). 

ii. / = on a convex sub-domain A(/) C f2: equation (MA) becomes degenerate 
on A(/) and T(f) — dA(f) is the free-boundary associated to this problem. 
The function / is C°° smooth on Q(f) = fi \ A(/) (cf. [CC]). The optimal 
regularity of / up to the interface will be discussed in this work. 

iii. / (-Pq) = at a single point Pq £ and f > 0, on fl \ {-Po}: equation becomes 
degenerate at the point Po- The function / is C°° smooth on \ {Pq} (cf. 
|CCj ). However, the regularity of / at Po is an open question. 



FULLY DEGENERATE MONGE AMPERE EQUATIONS 3 

We will restrict our attention from now on to the case (ii) above, where the 
solution / of (|MA|) vanishes on a domain A(/). 

By looking at radial solutions z = f(r) on fl = £?2(0) which vanish on i?i(0), we 
find that the expected behavior of / near the interface r — 1 is f(r) ~ (r — 1)1, with 
q given in terms of p by q — 53^- This motivates the introduction of the pressure 
function 

21 Q 
(1-1) 5 = 9 3 /', 



2-p 

A direct calculation shows that g satisfies the problem 



(MAP) 



with 



g det D-g + 9 (g 2 y g xx - 2g x g y g xy + g 2 x g yy ) = X{ s >o } in ft 
g = (p on<9f2 

n 1+ P 



2-p 

2 JL 

and (p(x) — q 3 93 1 . One observes that 9 > iff p < 2 which explains our assumption 
on p. 

A similar concept of pressure plays an important role in obtaining the optimal 
regularity of solutions to another degenerate equation, this time parabolic, the 
porous medium equation, namely the flow of a density function / of a gas through 
a porous medium given by 

(PME) f t = A/ m on M". 

The corresponding pressure g — / m_1 of the gas satisfies 

(1.2) g t = (m- l)gAg + \\7g\ 2 , on R". 

The pressure g is more natural in terms of the regularity. For a classical solution, 
the expanding speed of the free boundary dft(g) = d{g > 0} is |V<?|. If we observe 
that the free boundary expands with finite non-degenerate speed, g grows linearly 
away from the free boundary dfl(g), while the density / grows like a Holder func- 
tion whose Holder coeffcient depends on m, |CVW] . The pressure g is a kind of 



normalization of /. Then, g is Lipschitz on 1" [CVW], |CW) and smooth on fl(g) 

[Dm], m- 

A pressure-like function g = ^/2f for the parabolic Monge Ampre equation 
(1.3) f t = det D 2 f 
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has also been shown to be Lipschitz globally and smooth on f2(/) in |DH2] . [DL1J 
and [DL3] . 

Let us now turn our attention back to equation (|MAp . Our objective in this work 
is to establish the existence of a classical solution / of the problem (|MA[) . when the 
boundary data tp is such that the solution / vanishes on a region A(/) C f2 and 
therefore the equation becomes degenerate near the interface T(f) = dA(f). The 
concept of a classical solution will be discussed in section 12.11 To guarantee that 
such vanishing region exists, we assume that there is a classical super-solution ip of 
(|MA[) vanishing on a non-empty domain A(ip) C f2. In section [2~3l we will actually 
present an example which shows that this is indeed possible. 

Theorem 1.1. Assume that £1 C M 2 is a strictly convex bounded smooth domain 
and let tp G C 2 (fl), tp > on dfl and < p < 2. Assume that there is a classical 
super- solution ip of (IMAj) vanishing on a non-empty domain A{ip) C f2. Then, 
there is a classical solution of (|MAp and its pressure g, given in terms of f by 
(jl.ip . is C°° smooth on £l(g) up to the interface T. Consequently, f enjoys the 
optimal regularity f £ C k ' a with k = [gz^L a = — ^> ( C k ~ 1:1 , if k :— is 
an integer) and the interface T(g) is C°° smooth. 

A brief outline of the paper is as follows: in section [2] the concept of classical 
solutions of (|MA|) is introduced and the proof of its existence via the method of 
continuity is outlined. Section [3] will be devoted to the derivation of sharp a'priori 
derivative estimates for classical solutions of equation (|MAPp . These estimates 
play crucial role in establishing the C 2 s a regularity of classical solutions of (jMAPp 
which will be shown in section |4] (see in section l2~2l for the definition of this space). 
Based on the estimates in section [H we will conclude, in section [SJ the proof of the 
existence of a C^" up to the interface solution g of (|MAPp , via the method of 
continuity. We will also show that the pressure g is C°° smooth up to the interface. 

Notation: 

• fl C M 2 denotes a strictly convex bounded smooth domain in IR 2 . 

• p denotes a smooth strictly positive function defined on dQ. 

• For any g > on f2, we denote 

n(g) ={x€Sl\ g{x) > 0}, A(g) = {xeQ\ g(x) = 0} 

and T(g) = dA(g). 

• ds 2 denotes the singular metric defined in section 12.21 
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• llffllcgn = maxan(| Ajfll + |As| + .?)■ 

• C 2 ' a (fl) will be defined in section EJ and C^ 2+a {Q) will be denned in 
section [5] 

• v, T denote the outward normal and tangential directions to the level sets 
of a function g. 

• 9vi 9t) Qwi 9uti 9tt denote the derivatives of g with respect to v, r. 

2. Classical solutions and the method of continuity 

In this section we will define the concept of a classical solution of equation (|MA[) 
[ resp. of (|MAPp ] and sketch the proof of its existence via a method of continuity. 

2.1. The concept of classical solutions and the comparison principle. We 

consider the following generalization of equation (|MAP[) . namely 

(MAPh) g det D 2 g + 6 {g 2 y g xx - 2 g x g y g xy + g 2 x g vy ) = h X{ ff >o } 

where h £ C 2 (£l) and satisfies the bounds 

(2.1) < A < h < A" 1 < oo 

for some constant A > 0. 

We recall the notation fl(g) = {x € fl \ g(x) > 0} and A(g) = SI \ fl(g). On the 
free-boundary T(g) :— dA(g), where g = 0, we then have, from (|MAPh|) . 

(91 9yy ~ 2 9x9y9xy + 9 2 y 9yy) = glg T r = gin = h 

where v and r are inward normal and tangential unit directions to T(g) respectively 
and where n — grr/gu denotes the curvature of 

More generally denote by v, t the outward normal and tangential directions to 
the level sets of the function g. 

Definition 2.1. We say that g E C 2 {n(g)) iff 

Vgg 

extend continuously up to fl(g) and are bounded on fl(g). 
Define the non-linear operator 

(2.2) V\g\ := g det D 2 g + 6(g 2 x g vv - 2g x g y g xy + g 2 y g yy ). 
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Definition 2.2. Assume that g £ C ' 1 ^) H C%(Sl(g)) and that f = [q~3 g \ } 



q = 2^ is convex in fl. The function g is called a classical super- solution (sub- 
solution) of equation (jMAPhp if 



(2.3) 



V[g\ <h(>h) inCl(g) 

V\g\=BglK <h(>h) on T{g). 



The function g is called a classical solution if it is both a classical super- solution 
and sub-solution. 

If g satisfies (jMAPhp . then the corresponding convex function / = (q~i g^j , 
with q = 23^j satisfies the equation 

(MAh) det D 2 f = hf p . 

Definition 2.3. A convex function f is called a classical super-solution, sub- 
solution, or solution of (|MAh[) if the correponding pressure g belongs to C°' 1 (r2) R 



Cs(fl(g)) and is a classical super- solution, sub-solution, or solution of (jMAPhp 
respectively. 

Lemma 2.4. Let g\ be a classical super-solution and g<i be a classical sub-solution 
of (jMAPhp such that 

92 < 9i on <9f2. Assuming that £^(92) ^G?i)> we have 

92 < 9i in 0. 

Proof. Choose e > sufficiently small so that g\ := (I + e)gi < gi on dil and g\ is 
a strict sub-solution of (jMAPhl) . We claim that g\ := (I + e) gi < g\ in f2. 
Indeed, let us assume that touches g\ at a point Pq. If Po £ then 

h(Po) > V[ gi ](P Q ) > V[g e 2 ](P ) > h(P ) 

which is a contradiction. 

Hence, we may assume that Po £ 90(<?f). Clearly <9f2(gf) will also touch dft(gi) 
at Po- Then, at Pq, we have (g^v < i.9i)v and K2 < Ki, where k\,k,2 denote the 
curvatures of dft(gi), 90(^1) respectively. Thus at Po 

(Si)tt = (9i) v K i > {92) v K 2 = (sDtt 

and then 

7j ^ {gi)l{giU > (g £ 2)l(g £ 2)r 



h _ <~ ^2/„ \ ^ /„e^2/ e^ _ (l + e)h 



which is a contradiction. This finishes the proof of our claim. 
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Since (1 + E)gi = g% < g\ for any small e > 0, letting e — > we conclude that 
9i{P)<9i{P). □ 

Theorem 2.5 (Comparison Principle for Classical Solutions). Let g\ be a classical 
super- solution and gi he a classical sub-solution of (|MAPh[) such that 92 < 9i on 
dQ. Assuming that ^(#2) C £l(gi), we have gi < g\ in fi. 

Proof. The function g\ = (1 + 5 £ ) (172 — s)+ is also a strict sub-solution of (12.31) such 
that gf < 9i on f° r a small J e depending e. For large e > 0, fi(<7§) C fi(ffi) 
and let e decrease to zero. If (72 ^ 51 in 17, there is a positive e > such that 
<7f touches (71 at a point Po € ^(ffl)- The same argument as in the lemma above 
shows that Po can not be a point in fi(<7|), since g\ is a strict sub-solution. Also 
Po ^ dQ(g2); otherwise df^gf) would touch dQ(gi), for e > 0, which leads to 
contradiction similarly as in the proof of the previous lemma. □ 

2.2. The linearized operator near the free-boundary and sharp a' priori 
estimates. In sub-section !2.3l we will outline the proof of the existence of a classical 
solution of problem (|MAPp via the method of continuity. Our approach relies on 
the observation that one can obtain sharp a priori estimates for classical solutions 
g of the degenerate equation (|MAPh|) if one scales the estimates according to the 
natural singular metric corresponding to problem . 

To illustrate this better, assume that g is a classical solution of equation (|MAPh[) 
and that Po G T(g) is a free-boundary point. We will show in section 3, that 
g satisfies the a'priori bounds (|2.10[) and (|2.11[) near the free-boundary, which in 
particular imply the bound 

c< iDgl^cr 1 

for some c > 0. We may assume, without loss of generality, that g x > 0, g y — at 
Po so that it is possible to solve the equation z = g(x, y) near Po with respect to x 
yielding to a map x = q(z, y) defined for all (z, y) sufficiently close to Qo = (0, yo). 
The function q satisfies the equation 

(2.4) ~ z6cbD \ +eq ' qm ' =-H(z,y) 

where H(z,y) := h(x,y), x — q(x,y). Based on the a'priori estimates, we will 
show in section 4 that the linearized operator of equation (|2 .4[) near a function q 
satisfying the bounds (|4.8[) and (|4.9|) is of the form 

(2.5) L(q) = z an q zz + 2^/zai 2 q zy + a 2 2 q yy + bq z + cq 
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with (otij) strictly positive and b > v > 0. 

To apply the method of continuity one needs to establish sharp a' priori estimates 
for linear degenerate equations of the form (I2.5[) . These estimates become optimal 
when scaled according to the singular metric 

(SM) ds 2 = — + dy 2 

which is the natural metric corresponding this problem. 

Denote by B n the box B v = { < z < rj 2 , \y — yo\ < i] } and for any two points 
Qi = (zi,y\) and Q2 = (22,2/2) in B v , by s the distance function 

(DF) s(Q 1 ,Q 2 ) = \^-V^\ + \y 1 -y 2 \ 

with respect to the singular metric ds 2 . Let Cf(B v ) be the space of all Holder 
continuous functions on B^ with respect to the distance function s. Suppose 
that the function q belongs to the class Cf{B n ) and has continuous derivatives 
q t ,q z ,q y , q zz , q zy , q yy in the interior of B v , and that 

(2-6) q, q z , q y , zq zz , Vzq zy , q yy € Cf(B v ) 

extend continuously up to the boundary, and the extensions are Holder continuous 
on B n of class Cf{B v ) as before. We denote by C 2+a (B v ) the Banach space of all 
such functions with norm 

ll'?llc s 2 +°(8 r ,) = \\f\\c«(B v ) + II A?llc ? (B„) + lkt||c ? (8„) 

+ \\ z( lzz\\cf (Br,) + 1 1 Q»v II Cf(Br,) + hyyWcf (B„) ■ 



Definition 2.6. We say that g £ C 2 ' a (il(g)) if g is of class C 2,a in the interior 
of fl(g) and its transformation q £ C 2 ' a (Br)) near any free-boundary point Pq. We 
denote by \\g\\ c 2, a the corresponding norm. 

The following result follows as an easy modification of Theorem 5.1 in [DH2J. 

Theorem [DH] (Schauder estimate). Assume that the coefficients of the operator 
L given by (|2.5[) belong to the class C"{B n ), for some n > 0, and (a^) is strictly 
positive. Then, for any r < n 

ll?ll C |'"(B r ) ^ C ( M\\c°(B„) + \\Hcf(B„) ) 
for all smooth functions q on B^ for which Lq = h. 
The following result was shown in |DL2| . 
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Theorem [DL] (Holder regularity). Assume that the coefficients of the operator 
L given by (|2.5j) are bounded measurable on B^, Tj > 0, with (a-y) strictly positive 
and b > v > 0. Set dfi — x^" 1 dxdy. Then, there exist a number < a < 1 so 
that, for any r < r//2 

\\q\\c f (B r )<C (\\q\\c°(B v ) + (J B h'd^^j 

for all smooth functions q on B^ for which Lq = h. 

Based on Theorem [DL] and the sharp a priori bounds Theorem l2.8[ the following 
a priori estimate will be shown in section 4. 

Theorem 2.7 (C s 2 ' Q -estimatc). A ssume that g E C 4 (fl(g)) is a classical solution 
of problem (jMAPj) . with < p < 2, and that B p (0) C A(g). Then, there exists a 
constant C — C(\\(p\\c^ n , 0, p ) > such that \\g\\c 2 - a (n(g)) — C ■ 

Based on Theorem [DH] and Theorem l2.71 a C^' a solution of the problem (|MAP[) 
will be constructed via the method of continuity. It follows from Theorem [DH] 
and an inductive argument that the pressure g is C°° smooth up to the interface 
r(g), which readily implies that the interface is smooth (c.f. section 4). 

2.3. Existence of solutions via the method of continuity. We will now out- 
line the basic steps of the proof of the existence of a classical solution of (|MAP[) via 
the method of continuity. The proofs of these steps will be given in the following 
sections. 

According to our assumption in Theorem ll.il there exists a super-solution tp of 
(|MAp . i.e., -0 satisfies 

det(DV) < V, in ^ 

which vanishes on a non-empty domain A(ip) C f2. We define 

K detW) 

ipp ~ 

Before we proceed with the outline of the method of continuity, let us give an 
example which shows that there exist boundary values <f> for which such a super- 
solution can be found. 

Example. Set 

MP) = a(\P\ 2 -p 2 ) q + , i = ^ 
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and pick a c\ > so that 



2A < det(D*MP)) < 1 



^(P)f 2 

for some A e (0, 1) in ft(V>i)- When the boundary data <p in (|MA|) is such that 

tp > ipx, on 90 

we can modify ipi to a convex function il>(P), keeping the decay rate to zero on 
<9ft(V0, so that ijj(P) = ip(P) on dtt and 

Hence, ip is the desired super-solution. 

Going back to the method of continuity, we consider the following boundary 
value problems depending on a parameter t € [0,1]: 



(MAt) 



det(D 2 /(^)) = ((1 -t)h + t)P in ft 
f = ip on dil. 



Set /i t := (1 —t)h + t and observe that 

A < h t < 1 

since /i satisfies (|2.7[) . Hence, h t satisfies condition (|2.1j) . Also, since h t > h, a 
classical solution f(P;t) of (|MAt[) is a sub-solution of 



(2.8) 



det(D 2 f(P)) = hfP in ft 
/ = ip on <9f2 



while the given ^>{x) is a super-solution of (|2.8j) . Hence, by the comparison lemma 
[2751 if {ip{P) = 0} C {/(P;t) = 0}, then /(P;i) < ip(P) in ft. We are going to 
carry out the method of continuity starting with /o = ip(x) at t = 0, keeping 

W(P) = 0} C {/(P; t) = 0}, for < t < 1 

so that f(P;t) has a non-empty vanishing region A(/(P;t)), for every i € [0,1]- 
This justifies our assumption (H-2) below. 

Assume that / is a classical solution of (|MAt|) (we will drop the index t on / 
for the rest of the section). Then, the corresponding pressure function <?, defined 



(MAPt) 
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in terms of / by (II. satisfies 

g det D 2 g + 6 {g 2 x g yy - 2g x g y g xy + g 2 y g xx ) = h t in 
g — Tp on cT2 

for Tp = q^ify . 

We make the following assumptions: 
(H-l) n C Bi(0). 

(H-2) / and g vanish on a non-empty sub-domain A(/) = A(g) CC and 

B p (0) = {x e M 2 : |x| < p} C A(/), for some p > 0. 
(H-3) / is strictly positive and strictly convex on f2(/) = {a; € f2| / > 0}. 



(H-4) The pressure g satisfies g g C 4 (il((7)), i.e., in particular it is C 4 -smooth up 
to 90(c/). 

To simplify the notation, we will set from now on 

llsllcin = max(|A 3 .9l + \D ig \+g). 

In the next section we will establish sharp a-priori bounds on the first and second 
derivatives of the pressure g up to the interface dfl, as stated in the sequel. 

Theorem 2.8 (C 2 -estimate). A ssume that g is a classical solution of equation 
()MAPh|) in Q with < p < 2 and h € C 2 (f2) satisfying (|2 . 1 [) . Assume in addition 
that g satisfies the assumptions (H-l)-(H-^). Define the matrix 

y/ggvr 9rr 

with v , r denoting the outer normal and tangent direction to the level sets of g 
respectively. Then, there exists c — c(\\g\\ C 2^,\\h\\ C 2,9, X, p) > 0, for which the 
bounds 

(2.10) c < \Dg\ < c- 1 
and 

(2.11) c|^| 2 < ^-6^- < c" 1 ^! 2 , V^O 



(2.9) M = M 



hold on Q(g). 

Combining Theorem 12.81 with the Holder regularity Theorem [DL] , we will show 
in section U the following a priori estimate. 



12 PANAGIOTA DASKALOPOULOS AND KI-AHM LEE 

Theorem 2.9 (Cp^-estimate). Under the same conditions as in Theorem \2.8\ 
there is a uniform < a < 1 and C — C(\\g\\ C 2^ 7 ||/i||c 2 , 0, A, p ) < oo, such that 

\\9\\c 2 s - a (u(g)) ^ C - 

In addition the curvature K,(g) of the free-boundary T(g) is of class C a . 

The above result shows that the coefficients of the matrix (12.91) are uniformly 
Holder. This will be combined in section [4] with the Schauder estimate, Theorem 
[DH], to obtain the following regularity of g. 

Theorem 2.10 (Higher regularity). Under the same conditions as in Theorem 
\2.8\ and the additional assumption that h £ C°°(f2), the solution g of (jMAPhl) is 
smooth on £l(g) up to the interface T(g) which means that for every positive integer, 
there exists C]~ = C(\\g\\ C 2^ , ||/i|| C fc+2, 9, A, p,k) < oo for which 

\\9\\c*+ 2 - a tn&)) - Ck 

and the curvature n(g) ofT(g) is C k,a . It follows that g is C 00 -smooth up to the 
interface T(g) and that the interface is smooth. 

To implement the method of continuity, we next set 

I ={te [0, 1]| (|MAPt|) has a classical solution satisfying (H-l)-(H-4)}. 

Clearly / is nonempty since by the assumption of Theorem 1 1.1 1 ib is a solution of 
(|MAPt[) for t = 0. The existence of classical solution of (jMAPhp is equivalent to 
that 1 £ I. The method of continuity relies on showing that the nonempty set / is 
both open and closed in [0, 1] in the relative topology, which means that / = [0, 1] 
and hence 1 £ J. 

The closedness of / easily follows from Theorems 12.81 and 12 .91 as shown next. 
Lemma 2.11. The set I is closed. 

Proof. Let C / be a sequence converging to to- Then, there is a sequence 
of solutions {gk} of (MAPt), t = tk, and their free-boundaries T(gk) which have 
uniform estimates depending only on the boundary data and the domain Q. First 
we can extract a converging subsequence of the free boundaries to Tq and, 

among them, extract converging subsequence g^. converging to a function g$. The 
non-degeneracy estimate in (|2.10[) implies that Tq = T(go) and the uniform C\ ' a - 
estimate in Theorem 12.91 implies that go is a solution of (|MAPt[) with t = to. Hence 

t el. □ 
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The openness of I will be proved in Section [5] through the stability in the pa- 
rameter t, Theorem 15. 1[ which is similar to Theorem 8.5 in |DH2j . 

The method of continuity then implies the following existence of classical solu- 
tions. 

Theorem 2.12 (Existence of a classical solution). Under the assumptions of The- 
orem [Ll[ there is a classical solution g of li MAP]) which satisfies the estimates in 
TheoremslEMlKM and [2lU 

The rest of the paper will be devoted to the proof of Theorems 12.81 - Theorem 
12.101 which, in particular, imply Theorem ll.il 

3. Optimal Estimates 

In this section we are going to prove the optimal a'priori estimates stated in 
Theorem 12.81 We will assume, throughout this section, that g E C 4 (£l(g)) is a 
classical solution of equation (|MAPh[) in ft with < p < 2 and h 6 C 2 (J7) sat- 
isfying (|2.1|) . In addition, we will assume that g satisfies the assumptions (H-l) - 
(H-4) introduced in section l2~3l We recall the notation Q(g) — {x \ g(x) > 0} and 
h\\c 2 oa = maxanGAjfll + |Asl +3). 

We will first establish an upper bound on the first order derivative \Dg\. 

Lemma 3.1. Under the assumptions of Theorem \2.8l we have 



Proof. We set M := r 2 \Dg\' 2 = (x 2 + y 2 ) {g 2 x + g 2 ). We will show that M attains its 
maximum at dfl. This readily implies the desired bound, since r 2 = x 2 + y 2 > p 2 on 
f2(<?). (Notice that we cannot bound |-D<?| 2 from above by the maximum principle, 
if h 7^ 1 , so we need to multiply by r 2 ) . 

Let Po be the maximum point of M on il(g). Assume first that Po E 0,(g). We 
may also assume, by rotating the coordinates, that 

(3.1) g y = and g x > at Po- 

Also, since M x = M y = at Pq, we have 



max|U 5 | < C( P> 9,\,\\g\\ C 2,\\h\\ C i). 



(3.2) g xx = ^, and g xy = ^ at P 



which combined with (13. ip and (jMAPhp gives that 
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and hence 

,q,n r 4 h + y 2 gg 2 x 

(3 ' 3) r+< h ,«^, atPo - 

Let 

(3-4) A = {a ij ) = {gg ij +eg i g j ) t 

denote the transpose of the matrix G = (Gij) = {g gij + Ogi.gj)- This is the 
second order derivative coefficient matrix of the linearization of equation (jMAPhl) . 
Differentiating equation (|MAPh[) to eliminate the third order derivatives on ay My- 
and using (|3.1[) . (|3.2p and (|3.3p . we find, after a direct calculation, that 

6 

(3.5) dij My = y 



i=0 



with 



and 



P = Mr 5 (M9r - xg) and fe 6 = W(Qx 2 + y 2 ) 



N < (7(0,11%), i = i,...,5. 

Since r > x, assuming that M > Q~ x maxQ g we conclude that D > at Pq. Since 
the leading order term in (|3.5j) . when M is sufficiently large, is (b§M 6 )/D and 
&6 > we conclude that either M < C(p,p, A, HffHcg j II^Hc 1 ) a * Po or a ij > 0. 
In the latter case Pq cannot be a maximum point, contradicting our assumption. 

Assume next that Po £ T(g) and that M > at Po. We may assume again that 
(|3.1j) holds at Po, i.e. y is a tangential direction to T(g). Hence, M y = 0, M x < 
and M yy < also hold at Po- In addition, since g — at Po, equation (jMAPhp and 
(|3.1j) imply that 0g 2 g yy = h at Pq. We conclude, after some direct calculations, 
that 

(3.6) = 0, g xx <-~2-, 9xy = -^2-, 9w=Qjp atp o 

and 

2(x 2 -2u 2 )o 2 2r 2 /i 2 

(3.7) = 2 r 2 5:c g xyy + -i -/ 7 ^ + < 0, at P . 

On the other hand, differentiating equation (|MAPh[) with respect to x and using 
(j3~Tj) and ([32]) we find that 

, (1 + 2^)^ g g (1 + 2^ 



FULLY DEGENERATE MONGE AMPERE EQUATIONS 15 

which implies that g xyy = — + 20) g xx j g 3 . Substituting in ()3.7|) gives that 
_ 2 ( fe2 + V 2 ) M2 2 (! + 26)rxh + r 3 h x 2r 6 h 2 

which is impossible, if we assume that M is sufficiently large, depending on the 
data. This finishes the proof. □ 

We will next provide a bound from below on \Dg\. 

Lemma 3.2. Under the assumptions of Theorem \2.8[ we have 

\Dg\ >c(p,d,\,ma*\Dg\,\\h\\ c i)>0, on Q(g). 

ail 

Proof. For q > we set 

M := {x 2 +y 2 )- q {xg x + yg y ) = r- 2q g r , (x,y) € Sl(g) 
with g r denoting the radial derivative of g. 

Claim: There exists an integer q > 1 which depends only on data (on H/iHc 1 an d 0) 
and such that M > c(p,9,max a fi \Dg\, \\h\\ci) > on Q(g). Since, from condition 
(H-2) we have r 2 — x 2 + y 2 > p 2 for any (x,y) £ Fl(g), the claim readily implies 
the desired bound from below on \Dg\. 

We will next prove the claim by the maximum principle. Let Pq = (xq, yo) be an 
interior minimum point of M in ^l(g). We may assume, by rotating the coordinates, 
that 

(3.8) y = and x > p > at P . 
Since M x = and M y = at Pq we have 

xg xx - (2q - 1) g x = and x g xy + g y = at P 

and hence 

to n \ ( 2q ~ ^>9x , 9 V n , D 

(3.9) fitec = and g xy = at P . 

x x 

Substituting the above to equation (|MAPh|) . using also Q3.8p . gives 

,„ in . x 2 {l + h)+ g 2 y [g - (2q + l)6xg x ] 

3.10 g yy = " : at P . 

xg x [(2q - l)g + 0xg x \ 

Let A = (aij) be the matrix defined in (|3 .4[) . Differentiating equation ()MAPh[) and 

()3.8|) - (|3.10|) we find, after several direct calculations, that 

T - M -V4 ^ 

l, .— iviij — 2u i=0 — 
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with D = (M9x 2q + (2q - \)g) > and 

|& 4 |<C7(p > fl,A,max|X>5|,||fc||cO 

oil 

and 

6 = -(2q - 1) x- 2q - 2 g [2{q - 2)x 2 {h + 1) + 2qgg 2 y - x 3 h x }. 

By choosing q > 1 sufficiently large (depending on Ij/iHc 1 ) so that 

2(q-2)(h+l)-xh x >0 

we can make bo < 0. We conclude from the above that L < unless M (Pq) > c > 0, 
for some constant c = c(p, 9, A, maxgji \Dg\, \\h\\ci). This shows that an interior 
minimum of M must satisfy minM > c(p, 9, A, max^n l-Dfli) > 0. 

Assume next that Pq € r(g) is a minimum point for M. We may assume this 
time that (|3.1[) holds at Po. Hence, 

n r _ -2qxyg x + r 2 (yg y y + xg X y) _ n 

(J.llj iWj, - r 2(g+l) ^ U 

and 

/o 1 N n r [(1 - 2g) x 2 + y 2 ] ffa; + r 2 (yg xy + gg^) 

V 3 - 1 ^ _ r 2(q+l) - U 

and also, by equation (jMAPhj) . 9g 2 g yy = h at Po- Substituting g yy — h/(9g 2 ) in 
()3.1ip and solving with respect to g xy gives 

(3-13) *» = -teg + — ■ 

Substituting this in Q3.12[) and solving with respect to g xx gives 
^_ y 2 h , [pg-l^-pg + l)^ 



(3.14) 



'x 2 g 2 r 2 x 



Here we have used that x > at Pq. This follows from assumption (I3.1[) . (H-2) and 
the convexity of A(p). 

We next differentiate equation (|MAPh[) with respect to x and use that g = 0, 
g y = 0, g yy = h/{9g 2 ) and ()3.13|) to conclude that 

(3.15) g yyy = J^, at P . 

Also, we differentiate equation (IMAPh[) with respect to y and use that g = 0, 

gy = 0, .g yy = h/{9g 2 ), (|3.13l) - ()3.15[) . to conclude that 

r „,^ ^ ^ (l + 20)(l-2g)h 4(l + 2%V ff3; 

(3.16) 5 , ra < — + gj-g + — . 
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We next differentiate M twice with respect to y and use (I3.13| - (|3.16p . g y = 0, 
g = and g yy — h/(0g%.). We obtain, after several direct calculations, that 

M = x g xyy _ 2gx[x 2 + (2g - l)y 2 ]g x 2h + yh y 
yy r 2q r 2( 9 +2) + 0r 2q gl 

at Pq. Substituting p. 161) and = Afr 2? a; _1 in the above gives, after several 
calculations, that 



M vv < bi M + 



bo_ 

M 2 



with 
and 

b 



2q[2q + + 6q)y 2 - 6r 2 } 
h - 0^ 



x 2 ([1 + 40- 2q(l + 29)} h + 0yh y + 0xh x ) 



Q2 r 6q 

Observe that since h > A > we may choose q sufficiently large, depending on 
to make bo < — 1. Since, M yy > at P and r > p (by our assumption (H- 
2)) we conclude that M > c(p, 9, A, maxgo \Dg\, \\h\\c^) > 0, finishing the proof. □ 

We will next establish sharp upper bounds on the second order derivatives of g. 
We begin by an upper bound on the rotationally invariant quantity 

2„ r>„ „ „ i 2 2 „ 



G '■— 9x9yy ~ ^9xg y gx y + g y gxx — g v <jtt 

where v, r denote the outer normal and tangential directions to the level sets of g 
respectively. Since the level sets of g are convex (because the function / is convex) 
we have G > 0. 

Lemma 3.3. Under the assumptions of Theorem \2.8l the quantity G = g 2 g T r 
satisfies 

maxG<C(e,p,\,maxG,\\h\\ C 3)- 
n{ g ) dn 

Proof. Set M := G + \Dg\ 2 . We will estimate M by the maximum principle. Since 
g is assumed to be in C (Q(g)), and hence gg%j = at T(g), it follows from equation 
(jMAPhj) that M = 0~ 1 h + \Dg\ 2 at T(g). Hence, we only need to control M in the 
interior of £l{g). Assuming that the maximum of M is attained at an interior point 
Po E £l(g), we will show that 

i 4 

(3.17) ay M y - = M2n ^ n2 , £ M M\ at P Q 

i=0 



M 2 (l+g 2 ) 

with A = (dij) given by (|3.4|> . A4 > c(0, p, A, max^o G, \\h\\ C 2) > and |j4j| < 
C{0, p, A, max^n G, \\h\\ C 2), for i = 0, ...3. Since ay My < at a maximum point, 
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this will imply the inequality M < C(9, p, A, maxgo G, ||/i||c 2 )> showing that max^c,) M < 
C(maxgsi M, 6, p, A, maxgn G, \\h\\c^), as desired. 

To prove (|3.17l) . we begin by noticing that since M is rotationally invariant we 
may assume that (|3.1I) holds at P , i.e., g x > 0, = and M = g x (g yy + 1) at Pq- 
Also, by a standard change of variables (see in the proof of Proposition 4.1 in [S]), 
we may also assume that g xy = at Pq. Using (|3.ip we compute that 

M x = gl Qxyy + 2 g x g xx (l + g yy ) = 0, M y = g 2 x g yyy = at P Q 
implying that 

(3.18) g xyy = 9xx{ 9yy)_^ _ at 

Differentiating equation (lMAPh| in y, using (|3.18[) and solving with respect to g xxy 
we obtain 

(3.19) g xxy = at P 

9 9yy 

since g y = g xy = at Po- Also, differentiating equation (jMAPhp in x, using 
(|3.18|) - (|3.19|) and solving with respect to g xxx we obtain 

/o on N ^xffx + 9xx [(26 ~ g yv )gl + 2g (1 + g yy )g xx ] 

{6.ZU) g xxx = at f . 

9 9x9yy 

We next differentiate equation (jMAPhp twice in y, multiply it by g x and subtract 
it from ciij Mij to eliminate fourth order derivatives, while use p,18p - p.20p to 
eliminate third order derivatives. After several direct calculations, using also that 
g y = = g xy = at Pq, we obtain that 

4 



with 



and 



13 13 M2(l+ 5 2) 



A 4 = 36(1 + A6)g x 



\Ai\ <C(e,p,X,\\h\[ C 2,\\g\\ cl ), i = 0,---,3. 

By the previous two Propositions, < c < g x < C < oo. Hence, A4 > 0, while 
Ai, i — 0, • • • ,3 bounded. This shows that at an interior maximum point, M < 
C(6, p, A, maxgsi G, \\h\\c2), hence finishing the proof of the Lemma. □ 
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We will now bound 

(3.21) Q:=max(gD^g + 9\D ig \ 2 ), 

7 I — p 

from above, where the maximum in (|3.21[) is taken over all directions 7. Note, that 
in terms of the function /, we have 

Q = max(q 1/3 / i ^ £ /77) ) « = ' 
7 2 — p 

In particular, since / is convex, Q > 0. 

Lemma 3.4. Under the assumptions of Theorem \2.8l we have 

(3.22) maxQ<C(0,p,\\\g\\ci n ,\\Ho). 

Proof. We begin by observing that since g £ C 4 (£l(g)), by Lemma \3.1\ the bound 

Q = \Dg\ 2 < C(6,p,\, Nlcg n , I^IIcO holds. 

Assume next that the maximum of Q is attained at an interior point P$ € $l(g) 
and at a direction 7, so that 

Q(P ) =gD ll9 + 6\D l9 \ 2 . 

Let v, t denote the outward normal and tangential directions to the level sets of g 
respectively. 

Claim. Either Q(P ) < C(6,p,X, \\g\\c* , \\h\\ C 2) or 7 = v. 

To prove the claim, we begin by expressing the maximum direction 7 as 7 — 
\\ v + A2 t, with X\ + A 2 = 1 so that 

(3.23) Q(P ) = g [Xf g vu + 2 AiA 2 g vr + X\ g TT \ + X 2 g 2 v . 
Next, we use the equation (jMAPhj) expressed in the form 

(g guv + Q gl) 9tt = g gl T + h 

and the bounds in Lemmas 13.11 - l3~3l to first conclude the bound 

Q(Po)<C(9,p,\,\\g\\c* n ,\\h\\cz) 

unless gg V v is sufficiently large at Pq. If in particular gg vv > 0g 2 at Po, we then 
conclude from (jMAPhp and Lemmas ET21 and I3~3l that 

gglr < 1gg vv g Tr < C(9,p,X, \\g\\c% n , \\ h \\c 2 ) 99™ 

showing the bound g VT < C(6, p, X, \\g\\ C 2^, ||/i||c 2 ) y/9uu- Using once more the 
bound 9tt < C(6, p, A, \\g\\c 2 , IHIc 2 )i )> we readily deduce from f|3.23[) that Q(Pq) 
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becomes maximum when A2 = 0, provided it is sufficiently large, depending only 
on lls'llcg^) 0> P- This proves the Claim. 

If Q(Pq) < C(6,p,X, \\g\\c 2 i Il^llc 2 ) then the proof of the Proposition is com- 
plete. Otherwise, from the previous claim we may assume that Q(Pq) = g g vv + 9 g„ 
and also, since v is the maximum direction, that g g VT + 9 g v g T = at Pq, implying 
that g VT — at Pq, since g > and g T — at Pq. Also, by rotating the coordinates, 
we may assume that (|3.1[) holds at Pq, i.e., the direction of the vector v is that of 
the a;-axis. 

We will show that 

(3.24) > Oij Qij = g[(l + h)Q + C(9,p,X, \\g\\c% n , \\h\\c*)] at P Q 

with A = (a,ij ) given by (|3.4[) . Since h > 0, this implies the bound Q < C(9, ||/i||c 2 > \\g\ 
at Pq, which combined with Lemma T3 . 2 1 implies the desired estimate. 
To prove (|3.24p let us first summarize that 

(3.25) g x > 0, g y = 0, g xy = at P . 

Also, since Q = g g xx + 9g^. at Pq, equation (IMAPhp together with conditions (|3.25[) 
imply that 

(3.26) gyy^hQ- 1 at P Q . 

We next differentiate Q is x and y and use (|3.25[) to deduce the equalities 

Qx = g g X xx + (1 + 2 6) g x g xx , G y = g g xxy = at P 

from which we conclude that 

(1 + 2 9) g x g xx 

[o.z() g xxx — , gxxy — v, at _ro. 

Also, differentiating equation (jMAPhl) in x, using (|3 . 25[) and (|3.27[) gives 
(3.28) g xyy = ^ at P . 

We next differentiate twice the equation (jMAPhp in x to eliminate the fourth order 
derivatives from ay Qij and use p.27[l - (I3.28P to eliminate third order derivatives 
and also (|3.25l) - (13.261) to finally conclude, after several calculations, that 

(Hi Qij = g [(l + h) G + C(9, P , A, \\g\\c* a , \\h\\c) 

by the previous Lemmas and our assumptions. This finishes the proof of the Lemma. 

□ 
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We are now going to combine the estimates in Lemmas I3.HI3.4I to give the proof 
of Theorem El 



Proof of Theorem \2.8[ We begin by expressing (jMAPhj) in the form 

det M = ( g g vv + 9 sit ) g TT - g g 2 VT = h. 
Hence, it is enough to establish the bounds 

c < g guy + gt < c 1 and c < g TT < c" 1 

for a constant c = c(#, p, A, H^Hcrg j ||^||c 2 ) > 0. 

The bounds from above readily follow from Lemmas 13.31 and 13.41 combined with 
Lemma 13.21 The bounds from below follow from ( g g vv + 0g 2 ) g TT = h + g g 2 T > 
A > (from our assumption on h) and the corresponding bounds from above. □ 

We next re-state Theorem 12.81 in terms of the solution / of (|MAh| . 

Corollary 3.5. Assume that f is a non-negative weakly convex classical solution 
f of the boundary value problem (|MAh[) in f2, with < p < 2, which satisfies 
assumptions (H-1)-(H~4). Define the matrix 



(3.29) M - fa) 



1 l — 2p _ p 

q 3 f 3 Uu f 2 fur 
— e — i _ IlLe 

f 2 fur Q 3 / 3 frr. 



with v, t denoting the outer normal and tangent direction to the level sets of f re- 
spectively and q — 3/(2— p). Then, there exist a constant c — c(8, p, A, ||/ 9 ||c 2 n 7 ll^-llc 2 
for which 

(3.30) cl^/^^c" 1 ^ 2 , V^O. 

We will finish this section with the following lower bound on ^fg det D 2 g, which 
will be used in the next section. 

Proposition 3.6. Under the assumptions of Theorem \2.8l there exists a constant 
C = C , (||5|| C 2^, ||/i|| C 2, 6, A, p) > 0, for which the bound 

(3.31) yfg det D 2 g > -C 
holds on 0,(g). 

Proof. Set Z := (x 2 + y 2 ) ^fg det D 2 g. We will use the maximum principle to 
establish the bound Z > —C, which readily implies (|3.31[) . since x 2 + y 2 > p 2 on 
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Clearly, Z > — C on dfl(g). Assume that the maximum of Z is attained at an 
interior point Po G Since Z is rotationally invariant, we may assume, without 

loss of generality, that 

(3.32) g x > 0, g y = and g xy = 0, at P . 

Differentiating equation once and twice in x, y and using t licit Z x — Zy — at 
Po, we find, after several direct calculations, that at the minimum point Po where 
(|3.33p holds, we have 

1 3 

(3.33) < aij Zij = — g V At Z\ at P 

<ir A gg{g 22 

with A — (aij) given by Q3.4p . and 

A 1 = 13(x 2 + y 2 ) 2 g 4 1 g 22 + C^ 

and 

A 2 = -y/g (117 (x 2 + y 2 )^ g fg 22 + Cg) and A 



W„2„„ n .rin\ and A, = 



(x 2 + y 2 ) 2 g 22 ' 

The constants C = C(gi, g 22 , x, y) depend only on gi, g 22 , x,y and hence they are 
bounded, by Theorem 12.81 

We will show that a^ Zij < at Po provided that Z < is sufficiently large in 
absolute value and Po is sufficiently close to the free-boundary T(g), establishing a 
contradiction to a^ Z^ > at the minimum point Po of Z. 

It is clear from the estimates in Theorem 12.81 that A\ Z < and A 2 Z 2 < 0, 
provided Po is sufficiently close to the free-boundary 517, i.e. g is sufficiently close 
to zero. The term A 3 Z 3 is nonnegative, however we observe that 



A 3 Z 3 = - 7ty\* ' Z = CgZ 



Ax 2 g{^gZ) 2 
(x 2 + y 2 ) 2 g 22 



with C bounded, since y/g Z is bounded by the estimates in Theorem 12.81 Hence, 
Y^i=iAiZ l < A\ + A 3 < at P , provided that Z < is sufficiently large and 
Po is sufficiently close to the free-boundary T(g), which concludes the proof of the 
Proposition. □ 

4. C^ q -Regularity 



We will assume throughout this section that g £ C 4 (il(g)) is a classical solution 
of the boundary value problem (jMAPhj) in SI, with < p < 2 and h G C 2 (Sl) 
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satisfying (|2.ip . In addition, we assume that g satisfies the assumptions (H-l)-(H- 
4). Our goal is to establish a uniform estimate on the norm IMIc 2 .°(Q(gy)i as defined 
in section 2.2, by combining the a-priori estimates in Theorem 12.81 with the Holder 
Regularity result Theorem [DL]. We will obtain estimates which depend only on 
the data |Mlc§ n > \\Hc 2 ,S, X, p. 

Since the regularity theorem [DL] concerns with solutions on a fixed domain, 
we will first perform a change of coordinates, near the interface, which transforms 
the free-boundary problem (|MAPh[) to a nonlinear degenerate problem with fixed- 
boundary. The same coordinate change was used in |DH2| . We refer the reader to 
that paper for the detailed computations. 

Let Pq = (xq, yo) £ r(<7) be a free-boundary point. We may assume, by rotating 
the coordinates, that at the point Pq, 

(4.D "—isr*- 

Then, by Theorem 12.81 g x (P) > 0, for all points P — (x,y) sufficiently close to Po- 
Hence, we can solve around the point Pq, the equation z = g(x, y) with respect to 
x, yielding to a map 

x = q(z,y) 

defined for all (z,y) sufficiently close to Q = (0, y ). Using the identities 

1 „ % „ 1 „ 



9x j 9y — i 9x 

Qz q z 



,,3 



and 



. q y 1 \ 1 / Q y q y 

9xy = — ~ \~~^2 1 ZZ ' 1 z v I ' 9yy = I ~~2 Qzz ~~ 2— q zy + q yy 



qz \ qz qz J q z \qt i 

which yield to 



and 



2 i 2 \ 

9xx9yy ~ 9 X y = ~~a \1z Z qyy ~ q z y) 
Qz 



2 _ 9 i 2 — 

9y9xx L 9x9y9xy t 9 x 9yy — 3 qyy 

Qz 



we find that q satisfies the equation 

U 2 ) -zdctD 2 q + 6q z q yy = _ R 

qz 

with 

(4.3) H(z, y) = h(x, y), x = q(z, y). 

In addition, q is a concave function, since g is convex. 
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Consider the non-linear operator 

-z detD 2 q + 9 q z q yy 



Lq 



qi 



The linearization L of L around a point q has the form 

■Z q yy qzz + 2z q zy q zy + {9 q z - z q zz ) q yy 



L M - „4 

(4.4) 



4 z detD 2 q — 3 9q z q vv _ 

H e — qz- 

q\ 

Let us denote by B v the box 

(4.5) B n = {K<z<n 2 , \y-y \<v} 

around Qo = (0,yo) and by C"(B V ), C 2 ' a {B v ) the spaces defined in section 2.2. 
Our goal in this section is to establish the following result: 



Theorem 4.1. Assume that g S C 4 (fl(g)) is a non-negative classical solution of 
the boundary value problem (jMAPhj) on £1, with < p < 2 and h G C 2 (fl) satisfying 
condition (|2.1| . In addition, assume that g satisfies the assumptions (H-l)-(H-J^). 
Then, there exist constants !< a < 1, C < oo and r\ > 0, depending only on 
the data ||s||cg n ! \\h\\c^,0, X, p, such that for any free-boundary point Pq — (xo,yo), 
satisfying condition (|4.1[) , the function x = q{z,y) satisfies the estimate 

IMI C J+«(B„) < C 



on B v = {0 < z <n 2 , \y~yo\<n}. 



Consider the matrix 
(4.6) A={a lJ ):=q- i 
and the coefficient 



qyy V z qzy 
^zqzy 9q z -zq zz 



^ b = 4zdetD 2 q-39q z q yy 

qz 

A direct consequence of Theorem 12.81 is the following a-priori bounds on A and 

b. 

Lemma 4.2. Under the assumptions of Theorem \4-l\ there exist a positive con- 
stants c = c(\\g\\c<2^, \\h\\c2, 0) A, p ) and rja, for which the bounds 

(4.8) 0< C |^| 2 <^^^<c- 1 |^| 2 , V^O 
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and 

(4.9) < c < b < c- 1 
hold on the box B v , provided rj < r]o- 

Proof. By direct calculation 

(4.10) detA= zdeW2q : 9qzq ™ =h 

<A 

and 

(4.11) trA = [{g 2 y g xx ~~ 2g x g y g xy + g 2 x g yy ) + (gg xx + 9 g 2 x )] . 

9x 

By ([2J]), A < Act A < A -1 . The bound c < tr„4 < cr 1 follows from Theorem El 
and (|4.1ip . These two bounds yield to (g^J). 
Next, we observe that 

b = Az detD»g - 3 9 q z q yy = fe + ^ g2 
9z 

Theorem 12.81 shows that 6 < c _1 on B,. The bound from below 6 > c > on 
with 77 sufficiently small, readily follows from (|2.ip and (I3.31[) . □ 

We are now in position to show the uniform Holder bounds of the first order 
derivatives h y and h z of h on B, r 

Lemma 4.3. Under the assumptions of Theorem \4- 1\ there exists a number a € 
(0,1), and positive constants!] andC, depending only on the data ||<?||c? 2 n > ll^llc 2 ) ^> ^> 9 > 
such that 

\\qA\cf{Bv)<c and \\q y \\cf(Bv) <c 



Proof. We will first establish the bound for q — q y . Differentiating equation 
with respect to y we find that q = q y satisfies the equation L q (q) = H with 
H = —d y H, with L q given by f|4.4[) . Since d y H = h y + h x q y , using the notation 

H y (z,y) =h y (x,y) and H z (z,y) = h x (x,y), x = q(x,y) 

we conclude that q satisfies the equation 

(4.12) z«n q zz + ^Vzci 12 q zy + a 2 2 q yy + bq z + cq = -H y 

with onj and b given by (|4.6I) and (|4.7I) respectively and c = h x (x,y) = H z (z,y). 



In addition, Lemma 14.21 and our conditions on the function h, imply that equa- 
tion (|4.12|) satisfies all the assumptions of our C"-regularity result, Theorem [DL]. 
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Hence, there exists a number a in < a < 1, such that the Holder norm ||q ; ||c c, (Br 1 ) 
is bounded in terms of ||^||c°(B^) an( i II -^y 1 1 C°(B,,)- Since ||?||c°(B^) is uniformly 
bounded, the bound ||?j / ||c' t «(Ba) — ^ readily follows from our assumptions on the 
function h. 

We will now establish the C" bound for q = q z . Differentiating equation (j4.2|) 
with respect to z we find that q = q z satisfies the equation 

z qyy qzz - 2z q zy q zy + (z q zz - q z ) q yy 

(4.13) 



4zdetD 2 q- (39 + l)q z qyy „ «£„ 
e q z H e- = ni 

q b z q b z 

with Hi = d z H = h x q z = H z q z . We wish to apply the regularity Theorem [DL] 
shown in [DL2 to control the Cf norm of q = q z ■ However, our a-priori bounds in 
Theorem l2.8l do not imply that the term q zy /q z is bounded, since the bounds (j4.8[) 
only control \fz h zy . 

To control the C" norm of h z , we will apply Theorems 3.6 and Theorem 3.7 in 
[DL1] on certain super-solutions and sub-solutions of equation (|4.13|) . 

We begin by noticing that since the term q\ y jq\ is nonnegative, (|4. 13[) implies 
that q = q z is a super-solution of equation 

z qyyq zz ~ Zz q zy q zy + (z q zz - 9 q z ) q yy 

(4.14) 



qt 



4zdetD 2 q- {W + \)qzq V y . ^ „ 

s q z < t±i. 

ql 

Let us denote by (ajj) the matrix in (|4.6p and by 

4zdetL» 2 g- (W+l)q z q yy 

(4.15) b x := 

Hz 

and set 

(4.16) Li(q) := za lx q zz + 2^/zai 2 q zy +a 2 2q yy + biq z . 

A similar argument to that used in the proof of (|4.9p shows that bi satisfies the 
bounds 

(4.17) c<bi<c~ 1 , onB v 

with c = e(||p|| C 2 n , \\h\\ C 2,e,X,p) > 0. 

Following very similar computations to those in the proof of Lemma 5.9 in in 
|DL1] . we conclude: 
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• 9 = 1z is a super-solution of equation 

L\{q) <Hi, on B ri 

with Hi = H z q. 

• There exists a number j3 > 1, depending only on the a priori bounds, for 
which if (/i z — m) > on B v , for some positive constant m, then (72 := 
(h z — m)P is a sub-solution of the equation 

Li(q 2 ) > H 2 . 

• There exists a number /3 > 1, depending only on the a priori bounds, so 
that (73 '■= hP z is a sub-solution of the equation 

• There exists a number /3 > 1, depending only on the a priori bounds, so 
that for any constant M, q± := [M 13 ~ hP z ) is a super-solution of the equation 

Lx{qi) < H A . 

It can be shown, as in the proof of Lemma 5.9 in [DLlj . that the functions Hi, 
i = 1, .., 4 satisfy the bounds 

\Wih^)<C{\\g\\ci a ,\\h\\c'AX,p). 

The Holder regularity of the function h — h z on B v follows by combining the 
above with the Harnack estimate, Theorem 3.6, and the local maximum principle, 
Theorem 3.7 in |DL2| . along the lines of the proof of Lemma 5.9 in [DLlj . This 
yields to the bound ||&||c«(B f ) < C(||g|| C 2 f; , ||/i|j C 2, 6», A, p ). □ 

We will next combine Lemmas 14.21 and [4. 31 with the classical regularity results for 
strictly elliptic linear and fully nonlinear equations, to obtain the C^ ,a regularity 
of the solution q on the box B v defined by (14. 5|) around the boundary point Qo — 
(0,2/o, to); where Lemma l4~2l holds. 

Let Q r — (r 2 ,y r ) be a point in B v , where the index r indicates that the z 
coordinate of Q r is of distance r 2 from the boundary z — 0. For < /i < 1, denote 
by the disk = { z 2 + y 2 < /i 2 }. Define the dilation q r of q on D^, namely 
the function 



28 



PANAGIOTA DASKALOPOULOS AND KI-AHM LEE 



A direct computation shows that the function q r satisfies the equation 

~5 dct D 2 q r + 9 qlql n , 
(4.18) , zyy = -H r 

with z = 1 + z and H r (z, y) = H(r 2 + r 2 z, y r + ry). 

When P = (z, y) £ D^, with < (X < 1, then z > 1 - /i 2 > 0. It follows by the 
bounds of Lemma [4.21 and the bound 0<A<i?<A _1 , that (|4. 18[) is uniformly 
elliptic on D^. Hence, by the known results on the regularity of solutions to strictly 
elliptic fully- nonlinear equations (see in jCCj ). one obtains uniform C°° bounds for 
q r on _D M , in terms of ||g ,r ||L oo (Z3 (10 )) f° r an y < fi < fj,o < 1. Notice that, in 
addition, H^Hl^d^) is uniformly bounded, since q z is bounded in B n . The above 
discussion leads to the following lemma: 

Lemma 4.4. For any < fiQ < 1, there exists a constant C(/xo) depending also on 
llfllloln) \\ h \\c 2 ,P, A and p, such that 

for all < fi < fiQ. 



One may now combine Lemma 14.41 with Lemma 14.31 along the lines of the proof 
of Lemma 6.8 in |DL1] to establish the C° regularity of z h zz and ^fz h zy , as stated 
next: 



Lemma 4.5. Under the assumptions of Theorem \4-l\ there exists a number a in 
< a < 1 and constants C, rj depending only on the data ||<?Hc? 2 n > ll^llc 2 ) Pi-^tPj 
such that for any two points Qi = (z\, y\) and Q 2 = (22, y 2 ) in Br, we have 

\ziq zz {Qi) - z 2 q zz (Q 2 )\ + Iv^ifc(Ql) - yfz2qzz(Q2)\ < Cs(Qi,Q 2 ) a . 

Finally, the Holder estimate for q yy can be derived from the Holder estimates of 
q z , q y and z q zz , ^fz q zy and the regularity of H . 

Lemma 4.6. Under the assumptions of Theorem \4-l\ there exists a number a £ 
(0,1) and constants C, rj depending only on the data ||<7||c 2 n > ||^||c 2 > Pi\P> such 
that for any two points Q\ — {z\,y\) and Q 2 = (22,2/2) in Br, we have 

\%y{Ql) ~ Qvv(Q2)\ < s(Qi,Q 2 ) a . 

Following an inductive argument as in Theorem 7.3 in |DH2j . we can show higher 
regularity, as stated next. 
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Theorem 4.7. Assume that g £ C 2 ' a is a solution of (jMAPhj) which also satisfies 
the assumptions of Theorem \4-l\ and the additional assumption that h £ C fe+2 (f2), 
there exist constants < a < 1, C < oo and r) > 0, depending only on the 
data \\g\\c 2 A\^ , \\c k + 2 iPiKp> such that for any free-boundary point Pq = {xo,yo), 
satisfying condition (|4.1|) . the function x = q(z,y) satisfies the estimate 

lkllc*+-(*y ^ c 

on B v = { < z < rj 2 , \y — y \ < n } for any positive integer k. 

We are now in position to give the proof of Theorem 12.101 

Proof of Theorem \2.1(A Let 77 denote the uniform constant in Theorem l4.7l Consider 
the sub-domains 

n* v (g)={xen(g)\d(x,T(g))>r)} and Q v (g) = {x G Q(g)\ d(x,T(g) < r,}. 
The estimate in Theorem 14. 71 implies the bound 

\\9\\c5+«(nM) ^ c (h\\cl n A\h\\c^P:\p)- 
It remains to show that g £ C°°(Cl*(g)). Indeed, on f2*(g) we have 
< SoM < det D 2 f = hf p < C(A,max^) 

for a positive constants Sq and C(A,maxgo if). Hence, / satisfies a Monge- Ampere 
equation as those considered in |CKN] . 

The bounds in Corollary 13.51 imply the upper bound on any second derivative 
fa on fl*(g), and the lower bound of fa follows from the balance of the second 
derivatives det(£> 2 /) ~ 1 on fi*. Therefore / satisfies a uniformly elliptic equation 
and det 1/2 (D 2 f) is a concave operator. Hence, the C°° regularity of /, satisfying 
det D 2 f = hf p , on £l*(g) follows from the regularity theory for uniformly convex 
or concave fully-nonlinear operators ( |CC) ). □ 

The proof of Theorem \4-l\ readily follows from Lemmas 14.31 14.51 and 14.61 

5. Stability: / is open 

In this section, we will utilize the estimates of previous sections to show the 
following stability of solutions of (|MAPtl) in the parameter t. This will conclude 
the proof of the Theorem ll.il as discussed in section [231 



30 



PANAGIOTA DASKALOPOULOS AND KI-AHM LEE 



Theorem 5.1. Assume that go is a classical solution of (jMAPtjl for t = to, satis- 
fying conditions (H-l)-(H-^) and such that 1 1 c?o 1 1 o 2 ,c * — ^(ll¥'llc 3 n >l'>^>P) > Then, 
there is a S > such that for any t with \t — to\ < S, the problem (jMAPtl) admits 
a C 2,a -solution g{-,t). 

We will use the corresponding elliptic argument to the parabolic one which was 
used in section 8 of |DH2j . Since the two arguments are quite similar, we will only 
outline the proofs, referring the reader to |DH2] (see also in |DHlj ) for the details. 

We pick a smooth surface S, sufficiently close to the /o = (g~ 2 / 3 go) q , such 
that its inner boundary dS lies on the z = plane and its outer boundary is dfl. 
Denoting by T> a ring 



we let S : V 



V = {(u, v) E M 2 : I < u 2 + v 2 < 2} 
be a smooth parameterization for the surface S which maps 



d in V = {(u,v) : u 2 + v 2 = 1} to d ln S = SD{z = 0} and d out V = {(u,v) 
2} to d out S — dQ.. We can find a smooth vector vector field 



T 



(TA 

w 



which is transverse to the surface S n {z > 6} while it is pararell to the z — plane 
when < z < 5. Now we define the change of coordinate ip : T> — > R 3 by 



fx) 




(u\ 


y 


= f 


V 


V) 








Via this coordinate change, the solution z — f(x,y;t) of (MAt) will be mapped 
onto the graph 

-/ u \ 

V 

^ \w(u,v;t)J 

if z = f(x,y;t) is close to the surface S. By the choice of the parameterization S 
of S, we have 

(u, v) G d m V iff z = 0. 

In the other words, the interfaces T(g(x,y;t)) = d{(x,y) : g(x,y;t) > 0} will be 
always mapped to the fixed boundary d m V. 
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Definition 5.2. We say g(x,y; t) is of class C* ,2+a if the function w(x, y; t) belongs 
to the class C^ ,2+a (T>) . Finally, we say that g(x, y; t) are smooth up to the interface 
T(g(x,y;t)) if w(u,v;t) is smooth onT>. 

In addition, the equation (jMAPtp will be transformed to the boundary value 
problems 

!Mw(u,v;t) = (u,v)eV 
w{u, v] t) = tp(u, v) (u, v) € d out V 

where tp( u i v ) is the function, uniquely determined by (p(x,y), after the change 
of variables and Mw = F(D 2 w, Dw,w,u,v;t) is a fully nonlinear equation whose 
linearized equation at t = has the form (|4.16[) satisfying (j4.8[) . (|4.9[) . 

Theorem 15.11 follows by combining Theorem 8.4 in |DH2] and Theorems 14.11 and 
1471 
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